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A Model for Band-Gap Shrinkage in Semiconductors with 
Application to Silicon 


ABSTRACT 


A single mechanism is proposed as the dominant origin of energy-gap 
narrowing In heavily doped emitter and back-surface-field regions. The 
mechamism, Coulomb forces resulting In the correlated motion of majority 
carriers and screening by them, proposed before (by Sah), is here developed in 
more detail. The aim is l:o aid the understanding of energy-gap narrowing, and 
to provide a single formula for both heavily doped regions and also for 
regions in which holes and electrons exist together in large numbers (highly 
excited regions). Energy-gap narrowing influences conversion efficiency 
mainly through its effects on the open-circuit voltage. 


I. INTRODUCTION AND EXECUTIVE SUMMARY 


This report describes technical findings of work supported by contract 

no. 956525 for the period, October 1, 1984 to December 31, 1984. Details of 

the findings appear in the Appendix. In this section, we briefly summarize 
these findings and indicate their practical Implications for solar-cell 
design. 

1.1 General Background 

One of the tasks for this contract is the identification of key 
parameters that describe heavily doped regions. This identification together 
with an understanding of the associated physical mechanisms enables the 
possibility of systematic design. From a performance standpoint, the heavy- 
doping effects reduce the open-circuit voltage severely, thus also degrading 
the fill factor. There is some influence also on the short-circuit current— 

for short wavelengths from heavily doped emitters and for longer wavelengths 

for a heavily doped back-surface-field region. 

For the past ten years, energy-gap narrowing, commonly known in the 
literature as band-gap narrowing, has been put forward as one of the key 
parameters. It is known that the physical origin of band-gap narrowing is 
complicated. The detailed many-body problem has received attention from 
solid-state theorists. Their final results differ because of different 
approximations made in treating the quantum mechanical many-body problem. 
Moreover, their results apply strictly at T = OK, where the states relating to 
the impurity atoms are pure in the sense that they are not perturbed by 
lattice vibrations. 

To further complicate the problem, many solid-state theorists assume a 
periodic distribution in space of the impurity atoms. This ignores the 
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thermal diffusion— a random-motion process— that contributes to the placement 
of the impurity atoms in the lattice. Thus it ignores a random component of 
the spatial distribution of these atoms. 

On the other hand, other theorists emphasize this random component, which 
leads to band-edge distortion. This falls Into two catagories. 

*band tails 
*impurity bands 

The band tails have significance at higher impurty concentrations, whereas the 
impurity bands play a role for concentrations of the order of 10 18 /cm 3 . 
Experimental evidence exists to support these claims. Some workers have 
suggested recently that delocalized band tails or band tails having states 
with relatively high mobility are largely the origin of the band-gap 
shrinkage. 

It is highly important to identify the mechanism underlying band-gap 
narrowing. If one does not know the mechanism or mechanisms largely 
responsible, efforts to design processing to yield high conversion efficiency 
are frustrated. 

1.2 Our Findings Reported Here (see Appendix for details) 

We deal mainly with impurity concentrations above the order of 
10 18 /cm 3 . For such concentrations, experiment shows that the quantum density 
of states follows nearly a standard quadratic dependence. The states in an 
impurity band at lower concentrations move, for higher concentrations, into 
the majority-carrier band. Band tails exist, but we have evidence to believe 
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that they Influence mostly transport of minority carriers and have only a 
slight Influence on the energy gap for concentrations of the order of 
10^/cm^. 

Our findings appear In detail In the Appendix, a manuscript written by 
P,T. Landsberg, A, Neugroschel, F.A. Llndholm arid C.T. Sah. A brief overview 
fol 1 ows . 

A single mechanism Is proposed as the dominant origin of energy-gap 
narrowing in heavily doped and highly excited regions of silicon devices (and 
also of devices made of other semiconductors). The mechanism involves the 
correlated motion of the electrons and holes resulting from Coulomb forces. 
The resultant (Debye) screening reduces the energy of the system, thereby 
reducing the energy gap. In the literature, Sah first proposed this as the 
one mechanism that could produce true energy gap narrowing— that is, energy 
gap changes not arising from band edge distortion (bandtails and impurity 

bands). Our work fills in the details of Sah's proposal and tries to correct 
earlier attempts by others to develop a simple model involving screening. 

In addition to the development, a survey of earlier experiments, and of 
new ones obtained by us, is made. The agreement between the model and 

experiments on n- and p-type silicon is good as far as transport measurements 
near T = 300K are concerned. It is good in the sense that the discrepancies 
between theory and experiment are no worse than the discrepancies between the 
experimental results of various authors. It also gives a good account of 
recent optical (photoluminesence) determinations of energy-gap shrinkage 
at 5K. 

The most general form of the result of energy-gap shrinkage appears in 

Eq. (21) of the Appendix. For the degenerate limit, the model gives, in 

gausslan units 
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(3) 


AEq b 126.6(m*/m) 1/Z (e/10)' <2 (n/ xO 18 ) 1/6 

In regard to transport derived measurements of the energy-gap shrinkage, the 
model agrees best with those interpretations that do not assume equal minority 
and majority carrier mobilities. Elaboration on this remark appeared in the 
preceding Quarterly Report to JPL. 
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ABSTRACT 

A model for bandgap shrinkage in semiconductors is developed and applied 
to silicon. A survey of earlier experiments, and of new ones obtained by the 
authors, give an agreement between the model and experiments on n- and p-type 
silicon which is good as far as transport measurements in ~3O0 K range by 
various authors are concerned. In fact, the discrepancies between theory and 
experiment are no worse than the discrepancies between the experimental 
results of various authors. It also gives a good account of recent optical 
determinations of band-gap shrinkage at 5 K. 

The model is based on Debye screening and, apart from effective masses, 
temperature, dielectric constant and carrier concentrations, which must enter 
any theory of gap shrinkage in some form, it is parameter-free . The model 
gives in gaussian units AEg = 126.6(m*M} 1//2 (c/10)" 3//2 (n/10 18 )^ 5 meV in the 
degenerate limit. 


1. INTRODUCTION 


There has been much discussion of the problem of band-yap shrinkage in 
the recent literature. The problem is of considerable importance because the 
band-gap is one of the key parameters in semiconductor devices. Heavy doping, 
which leads to the effect, has been used increasingly as part of the drive 
towards microminiaturization. However, key aspects of the subject are still 
controversial. "Discrepancies in the values extracted from experiments" have 
been noted [1]. At the theoretical level there are also serious discrepancies 
which arise from the intrinsic complexities of many-body calculations and the 
difficulty one has in estimating errors. 'Some contributions to the band-gap 
shrinkage are not agreed upon as regards sign, while other contributions are 
estimated by some authors as large while others neglect them. We refer to 
important recent papers [1-3] for more details and for developments of the 
theory. Some time will clearly elapse before these matters are clarified, and 
this justifies a search for other approaches which are less ambitious but may 
capture some aspect of the essence of the problem, and lead to an adequate 
estimate. We pursue one such line here. 

Three points are to be emphasized at the outset. 

(i ) Simple Useful Formula 

Our main result for the band gap shrinkage AEq is equation (1), below. 
We observe that this formula represents the experimental facts to good 
approximation both at 5K (see Fig. 4) and at 300K (see Fig. 3). Thus the 
formula given can at the very least be regarded as a semi -empi ri cal summary of 
experimental results. It may be expected to .cover the experiments for the 
whole temperature range 5K - 300K and beyond. Section 4 discusses the 

comparison with experiment. The formula (1) applies to heavily doped 


semiconductors and to highly excited regions of semiconductors such as occur 
either under strong Irradiation or in the transition region of p-n junction 
subjected to a Urge forward voltage. The formula thus meets the need of many 
experimentalists and device engineers who require a simple analytical 
expression which can be introduced into other considerations that are 
sensitive to gap shrinkage. The device work is normally conducted at room 
temperature (while the r.vny-body theories are worked out for T=0). 

(11) Jell 1 urn Model 

What is the theoretical status of equation (1)? For its derivation the 

charged ions of the solid both of the host atoms and of the impurity atoms are 

imagined uniformly smeared out within the volume of the solid to supply a 
uniform background so that both ho,e and electron gases are electrically 
neutral. Thi s "jell i urn" model of two gases of interacting electrons or holes 
is of course well-known. But also allowed to Le present are .neutral atoms 

which, on being ionized, can give rise to electron-hole pairs. The long-range 
part of the Coulomb interactions leads to complicated and correlated motions 
among the particles of the two gases (which one may consider separately). The 
n-particle wavefunctions for these systems are not known. However, it is 
reasonable to assume that these correlations leave the short-range part of the 
interactions more or less intact, thus leading to short-range or screened 
potentials acting between largely independent particles, (the quasi -particles) . 
The screening parameter (equation (21), below) takes account of this 

situation. In such a model, interactions with ions, band structure, location 
of impurities, inter-valley scattering, etc. are clearly neglected, liven in 
T=0 many-body theories, however, these are problems whose importance has not 
yet been ascertained with certainty. Thus inter-valley scattering was 
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emphasized In £1], but lately ignored In [4], while the magnitude of the 

effect of impuf ty scattering on band gap shrinkage Is also not certain 

[2,4], Other unsettled questions of current research which arise In these 
treatments include questions of linear versus non-linear screening £3] or 

ordered versus disordered impurities on a lattice £1-3]. By ignoring these 

questions here we merely follow through the jellium (or screening) model. 

Our procedure has great simplicity, and it is afso highly approximate, 
and so does not compete with the much more sophisticated many-body theories. 
These are of necessity limited to T“0K at the present time, and so do not 
yield as yet what is needed by experimentalists. It is crucial to observe 
that we do not ignore these recent advances, £1-3]. Instead, we seek a direct 
and intuitive route to a formula which, while approximate, is likely to be 
useful. In fact, our paper represents a challenge, in that equation (1) 

should be deducible from many-body theories by appropriate approximations. 

In our derivation we go back to the Oebye-HuckcT theory of .electrolytes , 
as expounded by W. J. Moore [5] and as used in a related context by Morin and 
Maita [6] in a standard paper. Some generalizations of this early work 
/ (section 2) enable one to give the simple derivation required (section 3). 

(iii) Earlier Related Work 

Equation (1) has appeared before. Inkson [7] obtains it as his component 
and in the review [8] It appears similarly as equation (3.13). Equa- 
tion (1) was also used by Sah and col 1 aborators [9] to discuss experimental 
results. On the other hand similarly simple formulae, but different in 
essential details, appear elsewhere [10-12]. In the present paper it is 
however, given a central status and it is obtained by a simple derivation for 
the first time. 
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Our formula yields a larger value for AE(> for degenerate material than 
does the work [12]. The factor is of order 1.33 and can be quite significant: 
A typical difference is 200 meV compared with 160 ineV. The difference in the 
coefficients is due to the following circumstance. Both here and in [12] tK 
gap narrowing is attributed to a drop in an energy, X say, in the presence of 
screening by carriers, as against the energy X in the absence of screening. 
The energy X used for this purpose in [12] is the -energy stored in the 
field. As it is not obvious that this is the best choice, we chose in this 
paper for X the energy needed to create and separate ari additional electron- 
hoi 2 pair without imparting kinetic energy to either particle. This seems to 
correspond rather directly to the energy gap concept. The results differ 
therefore because the models differ in an important point of detail. 

2. THE MURIN-MAITA FORMULA 

\e 'eminence with a discussion of the formula to be used for band-gap 
shrinkage. It is 



AEq = e^X/e , 

(1) 

where e is the numerical 

value of the electron charge, 

e is an appropriate 

dielectric constant of the 
• 

material and 



x " "2^1 C6] 

(2) 

or 

X = Xjl [9] , 

(3) 


where 

= 4ne^(n + p)/ekT . 


(4) 


The quantity (3) Is essentially the Debye screening parameter . In [9] the MKS 
(S.I.) system was used, which has here been changed to c.g.s. (Gaussian) 
system by replacing the dielectric constant relative to vacuum by e/4n. The 
discrepancy between (2) and (3) will now be explained, 

Suppose the electrostatic potential due to a charge e can be written 

4>(e,r) *» b(r)e , 

where b is typically l/«sr» but can have other values (see equation (10), 
below). The electrostatic energy due to another charge e' at the appropriate 
point in this field is clearly bee', or for e' * e, 

U(e,r) s e<fr(e,r) = b(r)e 2 . (5) 


This is all that is needed to obtain the energy of the charge e in the 
potential d(e,r). The connection between U and aEq is made in equation 
(14). 

In [6] the energy calculation is approached by estimating the work done 
in charging an ion reversibly in its appropriate position r up to its total 
charge e. This is taken to be 


V(e,r) 


0 „ 0 
/ qd$(q,r) = ■c-rj-y / $(q,r)d<f>(q,r) 
q=0 q=Q 


IX e.r)] 2 _ U(e,r) 
2D(r) 2 


( 6 ) 


It was overlooked that this answer is half that which one would expect since 
it is assumed that the potential b(r)e also starts from zero, thus lowering 
the energy V below what it should be. 
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In Moore's treatment 


b(r) = \/e (independent of r) , (7) 

and it is obtained by splitting the normal Coulomb potential off the screened 
Coulomb potential: 

•»* 

-1 axpf-Nr) . -1 - 2£ . (8) 

The first term is the potential at distance r.due to a charge e in a medium of 
dielectric constant e. The second term is the potential due to the other ions 
in the electrolyte or semiconductor. The energy (5) or (6) (depending on what 
one wants to calculate) then gives a reduction whicn takes into account the 
electrostatic interaction among the current carriers. In [6] the combination 

of equations (2), (5) and (7) was taken to lead to the band-gap shrinkage 
(1). This was done by simply citing reference [5] where the second term in 
equation (8) was inserted in a free energy. Morin and Maita interpreted this 
as a change in chemical potential and hence in energy gap, and so applied the 
result as a correction to the np-product. The job was done in nine lines of 
text. 

It is useful, in order to bring out the semiconductor physics involved 
more explicitly, to indicate first a generalization of the above argument. 
One need not expand the exponent as in (8). It is perfectly valid to replace 
(8) by the exact resul t 

(e/er) exp(~\r) = e/er - (e/sr ) Cl - exp(-Xr)] . (g) 
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The second term of (9) has the same interpretation as the second term in 
(8). Th'.s (7) is replaced by 

b(r) = [1 - exp(-\r)]/er . (10) 

The generalized and corrected [for the factor 1/2 in (2)] Morin-Maita band-gap 
shrinkage is, using (10) in (5), 

U(e,r) = (e 2 /cr)[l - exp(-\r) J . (11) 

It depends on the distance r from the given charge for which the band 
shrinkage calculation is made. 

The introduction of concepts from electrochemistry [13] is not helpful 
for the present purposes, as it tends to distract from the essentials of the 
argument . 


3. THE PRESENT MODEL 

Vie continue to focus our attention on formula (1), with a view to 
arriving at it by a direct path suitable for semiconductor work. The first 
step now is to create in the jellium an electron-hole pair which is in a bound 
state for a very short time. The distance, a_ say, between the maxima of their 
wave packets is of the order of an Angstrom, so that the effect of the smeared 
out electron and hole densities will not affect the energy, U say, to create 
the pair. The normal Coulomb potential acts between the particles for r > a_ 
and is cut off at r s a_. Such cut-offs are often needed for small r as the 
Coulomb potential diverges as r *► 0. Imagine now the hole to be fixed at its 
instantaneous position, r = 0 say, and the electron to be removed to infinity 
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against the Coulomb attraction starting at the cut-off distance. (Alterna- 
tively one may regard the origin of tne coordinate system to move with the 
hole.) As remarked in section 1, the jellium consists of quasi -particles with 
short-range, or screened, potentials acting between them. The screening 
parameter (of equation (3)) is approximated as a constant (it really depends 
on the wave vectors involved in the Coulombic collisions). With these 

approximations the total energy supplied to create the .pair, and to separate 
it, is (Fig. 1) 


E Q (n,p) = W + (e 2 /ea)';xp(-\a) . (12) 

This quantity is interpreted as the minimum energy between a hole in the 

valence and an electron in the conduction band, where n,p are the carrier 
concentrations. We thus exlude the Burstein-Hoss shift and use in the 

sense of £gap,2 °f Mahan [3], If the semiconductor is highly nondegenerate, 
then screening can be neglected and (12) yields 

E g (0,0) = w + (e 2 /ea) . ( 13 ) 

As already explained, W and a_ are to be approximated as concentration- 

independent. By subtraction, 

AE g h E g ( 0,0) - E G (n,p) = (e 2 /ea)[l - exp(-Xa)] O U(e,a)] . (14) 

Thus one does indeed arrive at expression (11), but it is now rather more 
firmly linked to semiconductor concepts, a link having been established in the 
form aE q = U(e,a). Using Debye or Thomas -Fermi screening gives the same 

result in the limit of extreme degeneracy, namely for n + material 
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1/3 4e rn n 


1/3 

n 


» 


1 .e. 


(15) 



n 1/6 , n 1/6 . 

[ — a y] X 0.0863 A" x O 0.08(-4 tt) A" 1 for Si] (16) 

10 id 10 ia 


with e = 11.7 and for m n = m, where m n is the densi ty-af-states effective mass 
for electrons. 

The approximation 


exp(-Aa) ~ 1 - \a ( 17 ) 

requires the constraint \a < z wnere z ~ 1/4 (see Appendix). Hence from (16) 

,,,,1/2 ln 18 1/6 ,,,18 1/6 
a < (|Aj) (~ ) 11-592 A [= 12.53z[i!L-] A for Si.] 

n 

?n -l 

For n ~ 10 cm m Si, one has \ = 0.17A and a £ 1.45 A if z ~ 1/4A. 
Adopting (17), one finally obtains equation (1) from equation (14). 

Numerically, (1) and (15) yield 

m r ./m 1/2 1/6 

AE„ = 186.6 [— .!> ] [_" ] meV ( 18) 

b (e/10) 10 Ld 

where n is in ciiT^. For silicon with e = 11.7 and m n = m, this gives 

20 -3 

AE g ~ 215 meV at n ~ 10 cm . In fact, one. can put, using e = 11.7 for 
silicon, 
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(19) 


AE q - 


21S(n,„/n,) 1/2 [-^] 


1/6 


meV 


As noted In section 1, this is 1.33 times larger than a similar formula given 
by Lanyon and Tuft [12] for the degenerate limit. 

Although the band-gap shrinkage has been determined in (14) in terms of _a 
and K, a more sophisticated theory is needed to estimate a_. The importance of 
_a_ arises from the fact that it determines the relative contributions of the 
two terms in (12) and (13). The beauty of the present treatment is that no 
commitment needs to be made as regards the numerical value of - . The reason 
is that we need only (14) from which a_ cancels if (17) holds, so that the 
actual value of a., which enters only (12) and (13), is not required. 

The ratio of the two terms in (13) is 


_ e 2 /ea „ r 

= — " L 


eaE G (0,0) 




(2U) 


A divergence in R occurs when a_ is so small that W * 0. So a_ must be large 

enough to keep R finite. For n ~ IQ 30 cm -3 the limiting value of ^a is given 

by a_i 1.2 A for silicon. On the other hand, it must by virtue of (17) be 
small enough for \a < z, i.e. R must lie above the minimum value determined 
by z/\ a n 1//6 . Thus one can draw a curve such that, for given n, all finite 
R-values lying above this curve are permitted and all those lying below the 

curve are forbidden. .This curve is shown in Fig. 2 for Si and for z = 1/4. It 

shows that the theory holds for a very wide range of possibilities as regards 
the relative contributions of the two terms in (13). For n ~ 10 2u ciir 3 _a_must 
be less than about 1.5 A. 

The fact that the two quoted values ofjMie well within trie interatomic 
distance in silicon (5.42 A) makes the model internally consistent in the 
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sense that one does not expect screening to be of importance for distances 
smaller than the interatomic distance# 

The radio of the two terms in (12) is 

R' 5 , R _ 

Hence if R can lie above a certain value for given electron concentration n, 
then the freedom for R' is even greater in that it caq assume any value above 
a smaller number. 

In the Appendix it is shown that the theory begins to fail unless the 
electron concentration lies below a critical value v/hich is 

, ln 20 -3 
n c ~ 4 x 10 cm 


in the case of Si . 

The general screening length must make allowance for electrons and holes 
and has many interconnections with other parts of physics [14]. One finds for 
a parabolic band f 

0 2 a o 2 1/2 
aE G = T (eFT ^ N c F -1/2^ y c^ + f . 21 ) 

where CF_ 1 ^ 2 (a) = (l//n)J {x" 1//2 /Cl + exp(x~a)]dx) and y c s (F - E c )/kT , = 
(E y - F)/kT and F is the Fermi level. It is this formula which will be 
adopted here for all concentrations of holes and electrons. 


4. COMPARISON WITH EXPERIMENT 

Energy-gap shrinkage aEq obtained from our previous transport 
measurements on Si :As n-type layers at ~ 340 K [15] are shown in Fig. 3 as 
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full circles. A new result, obtained here, for Si : 6 p-type layers from 

transport measurements is shown as a full triangle. This new result was 
obtained using a p + -n-p transistor with uniformly doped 
(N^ 3 6 x 10 cm } 0.2 pm thick emitter. The measurement method was the 
same as in [15], The data was obtained in the range 320-380 K, Assuming aEq 
to be fairly temperature-independent in this range, it seemed reasonable to 
associate it with 340 K. The curve is based on our model, equation (21), and 
gives good agreement for m n /m » m /m - l.t and a majority carrier 

concentration n or p > 5 * 10 cm . 

For lower concentrations our model underestimates aEq. This is to be 
expected since the impurity band or tail electrons contribute to the 
screening, so that our estimate of K in (21) is expected to be an 

underestimate. At the concentrations well above 10^ cnT^ the impurity band 
is absorbed in the conduction band and one can return to a parabolic band. 

This sudden change near the metal -insul ator transition has been noticed for 

many other physical parameters [16]. 

In much of the literature Boltzmann statistics has been used to infer A£ G 
even if the semiconductor was degenerate. These data have been recalculated 
so as to isolate band-gap shrinkage and remove the effect of degeneracy. 

The results are shown in Fig. 3 and compared with equation (21). The 
experimental points are rather dispersed, but the general agreement is quite 
good. One reason for the deviations may be that some of the transport data in 
Fig. 3 were obtained from regions with nonuniform doping and the values of aEq 
reported are for the average doping concentration which is generally different 
from aEq for the uniformly doped samples. f The reason for this is that 
non-uniform doping gives rise to quasi -el ectric fields which leads to 
additional drift. In some other cases errors may have been introduced in the 
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evaluation of aE g by the use of geometrical corrections to the one-dimensional 
model . 

» 

The agreement between the model and experiment approaches the agreement 
between different experiments, and so can be regarded as satlsfatory. It is 
marginally better for the high concentrations (> 5 x lo^ 9 cm" 3 ) which are of 
the greatest practical importance for the heavily-doped regions of the 
emitters in bipolar transistors, in solar cells and in other semiconductor 
devices . 

The experimental estimate of Aliy by other workers tends to be somewhat 
smaller than ours for concentrations larger than about 5 x IQ 19 cm -3 . The 
reasons for this difference were discussed in detail elsewhere [23,24], but 
they are briefly repeated here. The aE g comes from the measured minority- 
carrier current which is a function of a product p[exp(AE G /kT)], where p is a 
minority-carrier mobility. Now if one uses a conventional assumption that the 
minority and majority carrier mobilities are identical, then one obtains aE g 
as reported by others and shown in Fig. 3. In contrast our analysis [15] 
assumes only a temperature independence of p in the measurement range, but it 
does not assume the magnitude of p, which is not well known [3,23,24]. Thus 
our interpretation of the transport data is more physical. The agreement 

between our transport data for the degenerate concentration range and the 
theory is excellent. 

Wagner's recent 5 K photoluminescence and photol uminescence excitation 
data [2b] are shown in Fig. 4. N- and p-type material appear to lie on the 
same curve of bandgap shrinkage against concentration. This suggests that at 
this temperature the optical characteristics are insensitive to this 
difference. We have therefore compared experiment with three theoretical 
curves based on (21), using the non-committal symbol m* for electron and hole 
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Vv„ 




effective mass. The actual values for heavily-doped silicon are not well 
known [16,26], but those for intrinsic silicon at 6 K are m n /m = 1.06 and 
nip/m a 0.59 [26]. These values have been included in the figure, though the 
best value for agreement with experiment is m*/m » 0.45. 

Some 2.4 K photoluminescence experiments have recently been used to infer 
the bandgap shrinkage appropriate to 300 K [27,28]. This data, corrected for 
Fermi-Dirac statistics, is given in Fig. 5. For intrinsic silicon at 300 K 
the effective masses of electrons and holes are m*/m ~ 1.1 [9,15,29]. 
Although the effective masses for heavi 1 y -doped silicon at this temperature 
are again not well-known, the model is seen to give satisfactory agreement for 
m*/m = 1.1* The transport data (Fig. 3) also shows good agreement between 
theory and experiment at ~340 K for m*/m ~ 1.1. The optical absorption 

experiments [30-32] have not been included in Figs. 6 and 6 as it is suggested 
that they underestimate aEq [25]. 

Lastly we wish to point out that the band-gap shrinkage in doped material 
has important contributions from one band only: the screening parameter k 
depends then on one concentration. For high optical excitation, however, both 
electrons and holes contribute to the screening parameters, which is thus 
substantially increased (see equation (21)). This results in a larger band- 
gap shrinkage, if the effects of any impurity-bands or tail states are 
neglected. This is illustrated in Fig. 6 which predicts the order of magnitude 
of expected experimental results. Figure 6 also shows the theoretical aE g for 
n-type GaAs. This is notably smaller than for Si since ('fl n )GaAs^ iT1 n^Si ~ 
0.066/1.1. The relevant experimental data is based on optical absorption {for 
a review see [8]) and is therefore not expected to be reliable. It; is not 
shown in Fig. 6. 
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5. CONCLUSION 


A model of band-gap shrinkage in semi conductors has been developed. Its 
application to n- or p-type silicon gives good agreement with both transport 
and optical measurements at 5 K and in the range 320^380 K. The discrepancies 
between theory and experiment are on the whole no greater than the 
discrepancies Inherent in the experimental determinations. Our model would be 
expected to underestimate AEq in the non-degenerate regime. Here the 
impurity-band and tail states are not merged with the majority band, and their 
contribution to gap shrinkage is not included in the model. 

The model outlined has two attractive features: 

(i) As it does not compete with models based on many-body theory it 
cannot claim to be as fundamental as other models recently 
proposed. But the benefit arising from this is its great 
simplicity. It can be applied to over a range of temperatures and 
carrier concentrations with considerable ease. 

(ii) It is essentially parameter-free. The only variables which enter 

are those which must be involved in any theory of band-gap 
shrinkage: effective masses, temperature, carrier concentrations 

and dielectric constant. (The latter would enter a many-body 
treatment in a different manner as a sum over electronic 
transitions .) 
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APPENDIX: CONSTRAINTS ON THE MODEL 


In order not to Introduce the parameter a_ Into the band-gap shrinkage 
formula, it is necessary to satisfy (17). Now the error Introduced in putting 

e"^ a a 1 - \a 


has the following values 


Xcl 

0.1 

0.2 

0.3 

% error 

0.6 

2.4 

5.9 


Let us put the requirement 


\a < z (z ~ 0.25) 

on the model. We also have from (13) that 

E 6 (0,0) > e 2 /ea . 

It then follows that 


\ < z/a < Ze E G (o,o)/e 2 

Some algebra now shows, using the degenerate limit, that one requires 


n 1 / 6 < n^ 6 = 1 10 ^ (J ) 1/5 . 

Here I 0 is 13.6 eV and a Q = 0.529 A is the Bohr radius. The critical 
concentration is 


n = (e/10) a 
C (m n /m) J 


(4z) 6 E r 6 x 0.659 x 1 Q 2U cm ‘ 
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where Eq Is expressed in eV. Using e » 11.7 and Eg ° 1.1 eV, the constra' 
is for silicon 


4.8O x IQ 20 cm" 3 

(%/m ) 3 


n < n c “ 
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Lu. 


FIGURE CAPTIOUS 


Fixj. 1 Schematic diagram showing the decomposition of the band-gap energy 

into W and the work done against attraction. More carriers are 
assumed present for curve 1 than for curve 2. 

ig. 2 A plot of R (equation (20)) against electron concentration n for 

Si. The values of R below the curve are i forbidden in the present 

20 -*3 

approximation, which begins to fail near n ~„n =4 [x IQ cm . 

L 

Fig. 3 Gap shrinkage as inferred from transport measurements for n-type 
layers from various sources at a mean temperatyre of ~34U K. The 

curves are based on equation (21), the upper curve is for in*/m - 1.45 
[10], the lower curve is for in*/m = 1.10, and e s 11.7 (Si) has also 
been used. The horizontal axis is the majority carrier 

concentration . 

Fig. 4 Comparison of (21) with recent optical data at 5 K [21]. 

Fig. 5 Comparison of (21) with recent photoeini ssion data at 300 K [23,24]. 

The data from [23,24] were corrected for degeneracy, as described in 
the text. 

Fig. 6 Comparison of two theoretical curves for Si at 300 K. based on 

equation (21). The lower curve applies to doped n-type silicon with 

m n /m = 1.4(5 [10]. The upper curve applied to intrinsic excited 
silicon with n = p, m n /m = 1.45 and m p / m = 1.08. Also shown, for 
comparison, is a theoretical curve for n-type GaAs with m n /m = 0.066 
and e = 12.5. 
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